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GENERAL INSTABILITY OF STIFFENED CYLINDERS 
By Herbert Becker 


SUMMARY 


Theoretical buckling stresses are determined in explicit form for 
circular cylinders with, circumferential and axial stiffening. The 
loadings are axial compression, radial pressure, hydrostatic pressure, 
and torsion. Analyses were confined to moderate -length and long cylin- 
ders. The investigation was based upon the use of a form of Donnell’s 
equation derived by Taylor which is applicable to orthotropic cylinders. 
The derivation of this equation is presented in this report. 


INTRODUCTION 


General Instability in stiffened circular cylinders has been inves- 
tigated for several cases of loading and types of stiffening. Analyses 
were performed utilizing energy methods, the three differential equations 
for orthotropic cylinders of Flligge (ref. 1), and various forms of 
Donnell's equation (ref. 2) applicable to orthotropic cylinders. 

Except for the analyses by Dschou (ref. 3) and Taylor (ref. 4), 
these investigations have not yielded explicit solutions for buckling 
stress. Furthermore, the analyses for external-pressure loadings have 
been confined to ring- stiffened cylinders. In the present report, 
explicit expressions are presented for axial- compression, external- 
pressure, and torsional loadings on moderate-length and long cylinders 
with both axial and circumferential stiffening. 

The axial- load solution Is that developed by Taylor, who utilized 
the simplifications employed by Donnell, together with the effects of 
orthotroplcity, and derived an eighth-order partial differential equa- 
tion for axial- compress ion buckling. In this report, the equation of 
Taylor has been extended to include torsional and pressure loadings. 

In addition, Taylor's assumption of a zero value for Poisson's ratio 
has been retained. 


The analyses are confined to moderate-length cylinders, for which 
the boundary conditions influence the buckling stress, and to long 
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cylinders, for which there is no boundary influence. Details of bound- 
ary conditions for the different loadings are discussed in the pertinent 
sections below. 

The investigation is restricted to elastic buckling. Linear theory 
is used in all cases. Section properties of the frames, stiffeners, and 
sheet are nominal, as depicted in figure 1. The actual properties to 
be used in design are discussed in part VI of the Handbook of Structural- 
Stability (ref. 5). This handbook contains a critical review of the 
field of general instability. It presents comparisons of theory and test 
data which delineate the utility both of the theoretical results presented 
herein and of other theories which do not employ the orthotropic-shell 
theory. 

A summary of the results of the investigations appears In table 1. 
Short discussions of the analyses of each case investigated are also 
Included. The derivations of the general instability stresses are pre- 
sented in appendixes A to D. 

This investigation was conducted at New York University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 


SYMBOLS 


Af 

Af 

A s 

A s 

a 

b 


d 

E 

F 

G 

% 


distributed area of frame including effective sheet. In., A^yd 
area of frame including effective sheet, sq in. 

distributed area of stiffener including effective sheet, in., A s /b 

area of stiffener including effective sheet, sq in. 
constant 

stiffener Spacing, in. 
frame spacing, in. 

Young's modulus, psi 
Airy's stress function 
shear modulus, psi 

distributed bending moment of inertia of frame, If/d, cu in. 
distributed bending moment of inertia of stiffener, I s /b, 


cu in. 
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If bending moment of inertia of frame including effective sheet, in. 

I s bending moment of inertia of stiffener including effective sheet, 

in.^ 

J = Jf + J s 

Jf distributed torsional moment of inertia of frame, Jf^d, cu in. 

J s distributed torsional moment of inertia of stiffener, J B /b, 

cu in. 

Jf torsional moment of inertia of stiffener, in.^ 

J s torsional moment of inertia of frame, in.^ 

ky buckling coefficient for unstiffened circular cylinder under 

lateral pressure 

kp buckling coefficient for unstiffened circular cylinder under 

hydrostatic pressure 

L length of cylinder, in. 

M x = EI b 

dx 2 


Mxy = - 1 GJ d 2 ^ 
^ 2 Sx dy 

M y = Elf ^ 


Sy' 


m,n parameters of axial and circumferential buckle wave lengths 
m = mL/jr 


N x axial normal load on cylinder, cr x t s , Ib/in. 

N-jjy shear load on cylinder, Tt, lb/in. 

Ny circumferential normal load on cylinder, Oytf, Ib/in. 
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p pressure loading on cylinder, psi 

Qx,Qy transverse shears on shell element, lb 

R radius of cylinder, in. 

T torque loading on cylinder wall, in. 

t thickness of cylinder wall, in. 

t effective thickness of cylinder wall in shear-buckled Btate 

tj distributed area of frame, Ap/d, in. 

t s distributed area of stiffener, A g ^b, in. 

u,v,w displacements in x-, y-, and z-directions, in. 

Wjjjh deflection constant 

x axial coordinate, in. 

y circumferential coordinate, in. 

z radial coordinate, in. 

Z L parameter for unstiffened cylinder, L 2 /Rt when v = 0 

(3 = n/m 

7 shear strain 

e x axial normal strain 

e y tangential noimal strain 

?\ half-wave length of buckle in circumferential direction, in. 

v Poisson's ratio 

Pf radius of gyration of frame section, p^ 2 = Ifjkf = I fjt f , in. 

o c axial- compressive -buckling stress, psi 

<*x 


axial noimal stress, psi 
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T 


circumferential- compressive -buckling stress; also, general 
circumferential normal stress, psi 

shear-buckling stress; also, general shear stress, psi 


HISTORICAL BACKGROUND 


Early theoretical investigations into the buckling of stiffened 
circular shells were performed by FlUgge (ref. l), Dschou (ref. 3), and 
Taylor (ref. 4 ). Flilgge derived the trio of linear equilibrium equa- 
tions analogous to those for isotropic cylinders which have been utilized 
by many investigators. Dschou solved these equations for stiffened 
circular cylinders under axial load. 

Taylor derived a differential equation for axially loaded ortho- 
tropic circular cylinders utilizing the same approach as did Donnell 
in obtaining his well-known eighth-order partial differential for iso- 
tropic circular cylinders (ref. 2). The result is a relatively simple 
explicit solution for the general instability stress which reduces to 
the classical result for an isotropic cylinder. There Is a slight dis- 
crepancy due to Taylor's assumption of a zero value for Poisson's ratio. 

Buckling of ring-stiffened circular cylinders under external pres- 
sure was investigated by Salerno and Levine (ref. 6) and by Kendrick 
(ref. 7) by using the energy approach. Kendrick's calculations led to 
higher Instability stresses than those obtained by Salerno and Levine, 
who omitted certain terms in the energy equations. 

Bodner investigated this case (ref. 8) using an equation of the 
Donnell type similar to that presented in this report in which Poisson’s 
ratio is included. The complete length range was analyzed. 

Stein, Sanders, and Crate determined the section properties of ring 
stiffeners required to avoid general instability in cylinders loaded in 
torsion (ref. 9 ). Hayashi (ref. 10) analyzed stiffened cylinders in 
torsion employing the identical approach used by Donnell for isotropic 
cylinders. Hayashi showed that Donnell's results are a li miting case 
of the buckling of orthotropic cylinders in torsion. The complete length 
range is included in this treatment, and the results are obtained in the 
same form as found by Donnell, except for the more general form of the 
parameters . 
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THEORETICAL RESUITS 


The theoretical buckling stresses for orthotropic cylinders -under 
various loadings are shown, in table 1. Derivations of these expressions 
are presented in the appendices. Buckling stresses were determined on 
the assumption that the spacings of the longitudinal stiffeners and the 
circumferential frames were small enough to consider the cylinder to act 
as a uniform orthotropic shell. Effects of boundary conditions on the 
general instability stresses for the various loadings are discussed in 
the following paragraphs. 


Axial Load 

The solution chosen by Taylor for the case of axial loading repre- 
sents the waveform assumed in the classical solution of simply supported 
isotropic cylinders, in which linear theory is used. It is applicable 
to both moderate -length and long cylinders. 


External Pressure 

The solution for external pressure is identical to that for axial 
load and was applied by Batdorf to external radial and hydrostatic pres- 
sures on simply supported cylinders (ref. 11). For long isotropic cyl- 
inders Donnell's equation leads to a buckling stress that is too high, 
compared with the results of an exact analysis, by a factor of 4/3. 
Although the corresponding result is obtainable for orthotropic cylin- 
ders by the method presented in this report, Levy's result for radially 
loaded rings (ref. 12) has been presented instead Bince it corresponds 
to the exact solution. 

The analysis is simplified for external-pressure loadings by assuming 
that the ratio of axial to circumferential wave length of the buckle is 
negligibly small compared with unity. This simplification is justifiable 
for moderate -length cylinders. Furthermore, the explicit result agrees 
with Bodner’s data for ring-stiffened cylinders in this length range. 


Torsion 

The solutions and simplifications applied to stiffened cylinders 
loaded in torsion were employed by Gerard and Becker for moderate -length 
isotropic cylinders (ref. 13). The long-cyli n der solution does not sat- 
isfy boundary conditions, although it apparently is a satisfactory repre- 
sentation of the buckle waveform for isotropic cylinders. However, in 
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this length range a satisfactory solution to the problem may be obtained 
although boundary conditions are not satisfied for isotropic cylinders. 
The same situation was assumed to apply to stiffened cylinders. 

The moderate-length solution satisfies the requirement that w = 0 
on the boundaries but does not correspond to vanishing moments at these 
locations. However, the buckling stress for an isotropic moderate -length 
cylinder in torsion is relatively insensitive to specifications on these 
latter quantities, since there is less than 10-percent difference between 
the theoretical buckling stresses for simply supported and clamped edges. 

The general instability behavior of orthotropic cylinders is con- 
sidered to parallel the buckling behavior of moderate -length and long 
isotropic circular shells. Consequently, the same, solution to the buck- 
ling equation and the same simplifications in the mathematics are assumed 
to be applicable. The mathematics was simplified by assuming that the 
ratio of circumferential to axial wave length was negligible compared 
with unity. 


College of Engineering, 

New York University, New York, N. Y., 
April 10, 1956. 
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APPENDIX A 

DERIVATION OF DIFFERENTIAL EQUATION 


The derivation of Taylor's differential equation in which the 
effects of shear and circumferential normal stress are also included is 
presented below. The geometric properties of a stiffened cylinder are 
portrayed schematically in figure 1. 


Taylor derived the equilibrium, equation for an orthotropic cylinder 
by combining the compatibility equation for forces in the plane of a 
plate element with that for equilibrium of forces normal to a plate ele- 
ment with initial curvature in one direction. These two situations are 
depicted in figure 2. 


For equilibrium in the "plane" of the element 



( 1 ) 


and for equilibrium normal to the element 


a 2 ™* _ 2 Atjr ^ 


Sx2 


dx by $y2 


b^v 

N x — — + 
dx2 


2N- 


xy 


b^w 
dx by 



~ + P = 0 (2) 


in which i + is the sum of the curvatures in the y-direction due 

R By 2 

to the initial shape of the shell and to the deformation under load. 

It is reasonable to assume that Poisson's ratio vanishes for a 
stiffened cylinder fabricated in the usual manner of aircraft construc- 
tion, in which case Hooke's law and the stress-strain displacement 
relations assume the forms 
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= 

dx 


El = 5l_ = 

E Etf 


e v = 


bv 

¥ 



I = Em - y = du , dv 
G Gt by bx 


Mv = El. 


d^w 

ax 2 




M xy = 


1 

2 


GJ 


b^w 
bx by 


(3a) 

(3b) 

(3c) 

(3d) 

(3e) 

(3f) 


From equations (3a), (3b), and (3c), the compatibility relation 


&ge x , a S a a r 1 

dy 2 dx 2 H dx 2 

follows, which, is equivalent to 

i 1 _ _1_ ^^xy _ 1 3^w 

Et s 3 y 2 Etf ^ x 2 Gt dx by R 3 X 2 


00 


(5) 



10 


NACA TN 4257 


If the Airy stress function P is introduced 


N 

x dy 2 

N = 8% 
y 


N - d 2 P 
^ dx dy 


then equation ( 5 ) he comes 


1 d^F 1 d**P _1_ dV _ 1 dfw 

. ®bf di? 0^ dx^dy2 Ets dy^ - R dx^ 


( 6 ) 


(7) 


This is the compatibility equation for forces in the plane of the 
element. The equilibrium equation for forces normal to the plane of 
the element may be obtained from equation ( 2 ) by writing the term 


N = i^ +N XJL 

dy 2 Rdx2 ^-2 


d^w 

dy 2 


and utilizing equations (5&)> (5e),-and (5f) which yield 


EI e 


dV 


dV 


d^Sr 


d^ 


d^w 




hL" + GJ + Elf Ji-f + N x ^ + 2Njry v ■ < - + N v 

5? dx 2 dy 2 dx 2 ^ ^ dy * by 2 


+ P 


1 
R 


(8) 


The function F may easily be eliminated from equations ( 7 ) and (8), 
which will lead to the relation 


cv 10 ^ 
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e_ a 4 i_ a{{A 
Gl Sx 2 ^ 2 t s by k J 



a 4 + gj a 4 + 
ax 4 ® a x 2 ay 2 


if 



+ 



+ 2N- 


xy ax ay + 



1. Ar 
2 ax 4 


( 9 ) 


If it is assumed that p is constant. For an isotropic cylinder this 
will reduce to Donnell's equation for v = 0. 
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APPENDIX B 


GENERAL INSTABILITY UNDER AXIAL COMPRESSION 


Taylor utilized, as the solution to equation ( 9 ), 


w = w^ sin mx sin ny 


( 10 ) 


for the analysis of axial compressive buckling of an orthotropic circular 
cylinder. For this case 


N x = 0 x t 8 


N y = Nxy = 0 


( 11 ) 


Upon solving for a x and minimizing with respect to m and n. 


2E 

° c = Rt s 


p 2 I f + (GJ/2E) 


(p 2 /ts)+ (E/2G?) 


1/2 


( 12 ) 


and 


l/m = R 


= -R 1 / 2 


I B + P 2 (GJ/E) + p 4 I f (l/t f ) + (p 2 E/Gt) + (p^Aa) 


1/4 


( 15 ) 


where 
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where 


t 

t f I s “ t s I f 

2t f 

. 

t s i f - (g^j/e 2 ) 


Vs - 

'g^j/e 2 ) 

^ - 

[g^Ej/e 2 ) 


The half-wave lengths are it/m. and jt/n in the longitudinal and 
circumferential directions, respectively. 

When t s = tf * t = t and I B = If = t^/l2, then p 2 = 1 and 
equation (12) becomes 

= (3) _1 / 2 Et/R (15) 

which is the classical equation for buckling of a long circular iso- 
tropic cylinder under axial, load when Poisson 1 b ratio is zero. When 
v = 0.3, the error is 4.6 percent. 
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APPENDIX C 

GENERAL INSTABILITY UNDER EXTERNAL PRESSURE 


For external-pressure loading, equation (9) becomes 


1 a 


+ 4 


i ^ 


Sx^ ax 2 8y 2 t 6 kj 


T b k GJ b k 
J-s — r + 


+ I f + 

Sx^ E 8 x 2 8 y 2 Sy^ 


ax 2 


f N * ^ + » 


a 2 N 

y Sy 2 ; 


v + 1 ^ = 0 

r 2 a? 


(16) 


For uniform radial pressure, N x = Nxy = 0, and 


Thus 


N y = pR 


jl a 4 . e a 4 , i a» \ L tf_ + 
a? Gt a x 2ay 2 ■ts 8y4y \ B 8x4 


(17) 


a + i f St + a ^Vv) i aV . o 

E dx^dy 2 dy^ E dy 2 / P . 2 dx^' 


(18) 


Let 

w = v mri sin mx sin ny 
Then equation (l8) becomes 


Eka 2 n 2 n^'\/ x 4 GJ p p T 4 Rn 2 ) 

V*f Qt ^b)\ E f Ej r 2 


(19) 


( 20 ) 
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and consequently the pressure is found in the form 

+ 91 m 2n2 + If n 4 + ^ 

E E R 2 /s^ + Ec& 2r >- 2 ' + sJt 

\t f Gt t s 

or, with p = n/m, the expression for the pressure becomes 


+ (k^ + ; | p 4 + |£' 


tf G? 


Moderate-Length Cylinder Under Radial Pressure 

The critical pressure for a moderate- length cylinder under radial 
pressure is obtained by minimization of pR, since R is constant. 

From M = 0, 
dm 

/ »~l/4 / v -l/4 


i = (e) 1/2 (i 6 + f p 2 + ijA (± + *£<.£: 

m \ E /\tf Gf t s> 


which is the same as for an axially compressed orthotropic circular 
cylinder. When equation (23) is substituted into equation (22), it is 
found that 




(B|3) 2 1 1. ^ 2 £ 

\tf Gt H t B 


This is equal to a c yp 2 , a n d becomes a c when 3 = 1. 
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Actually m = — and n = rt/X where X is the half-wave length 
L 

in the circumferential direction. Then equation (22) becomes 


f ■ (f) 2 (| + f + vj + 



( 25 ) 


For t he min imum integral value of a in the solution of equa- 
tion (25), a = 1. Then with P » 1 (where (3 = L/X), which assumes 
that the axial wave length is much larger than the circumferential wave 

length. 


= 

E 



.-1 




(26) 


This is obtained by neglecting terms of lowest order in 

frcM $k£l.o, 

dp 


p. Then 



( 27 ) 


and thus 


E L\ 




Since 


cr, r = 



then 


a y 


5-51 



(28) 


( 29 ) 



MCA TN k237 


5 - 51 1 mm i± 


l \t f A t f y i r 2 


= 5 . 51 E 


. lA n 3/2 

M M 1 


For isotropic cylinders , 


, ,?A / 5 \ l/2 

5.5! l/iA [£) 

L 12 \B / 


Using the notation 


3/0 

'V = °- 85 E (t) 


ky«%Jt 2 

12(1 - v 2 )l 2 

t 2 , 

Then with. Z L = — , equation ( 31 ) reduces to 
Rt 

k y = 1.033Z L 1//2 (3; 

(note that v = 0) which, checks figure 1 of Batdorf *s report (ref. 11), 
in which. 


ky = 1.039Z L 1 ^ 2 


for Zj, > 100. This result can he obtained analytically. Fran Batdorf ' s 
report, for radial pressure on an isotropic cylinder. 
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k y = 


_U 


“ 2 ±j££ + 


12Z L 2 m l4 ' 


P 2 n 4p2^ffi2 + g2j< 


where m = mit/L. For hydrostatic pressure. 



When m = 1 and p 2 »1, 


^ - “p - <> 2 + ^ 


Minimization of equation (36) leads to equation (33)- 


(34) 


(35) 


(3 6) 


Moderate-Length Cylinder Under Hydrostatic Pressure 

For the moderate-length orthotropic cylinder under hydrostatic 
pressure, equation (25) becomes 


pR/m^ 
E \2 


+ 



I B m4 + m 2 n 2 + Ifn^ 
E 


+ 




2/m^ t Em^n 2 t n^\ 
rtf Gi t s i 


(37) 


or 


PR 

(p 2 + 

= m 2 (l s + — p 2 + Ifp^) + 

R2m 2 (i- 

+ 4 

p 2 

E 

\ 2 I 

\ E / 

\tf 

Grb 



(38) 


When p 2 » 1, the results of hydrostatic and radial, pressure are 
equal, and equation (30) will apply. 
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Long Cylinder Under Radial Pressure 

For a long cylinder under radial pressure, the buckling mode 
sponds to that of a ring, for which Levy obtained the result 


P = 



With 



it follows, from equations (39) and (40), that 

«V - 3 E N r ) 2 


corre- 

(39) 


(40) 


(41) 
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APPENDIX D 

BUCKLING OF CIRCULAR CYLINDERS UNDER TORSION 


Moderate -length. Cylinders 


For shear loading of moderate -length circular cylinders 


N^y = tr = T/2*R 2 
N x = Ny = 0 


(42) 


As was shown by Becker and Gerard for isotropic cylinders of moderate 
length, a useful solution for r cr is obtained when the expression 


w = sin *** + HZ - sin 
\ L R/ 


(m + 2)^ + 


2 £ 

R 


; ] 


( 43 ) 


is used in equation (9). This satisfies w = 0 at the cylinder ends, 

although dw/dx and d 2 w/dx 2 are not prescribed. Such freedom from 
boundary conditions is indicated for isotropic cylinders by the small 
difference (less than 10 percent) between the buckling stresses for 
clamped and simply supported cylinders of moderate length. 


Equation (43) is adapted here for orthotropic cylinders. When it 
is substituted into equation (9); with the additional stipulations that 

(mrtR/nL) 2 « 1 and j(m + 2)itR/nL] « 1, then it is found that 


LI f 

fnf Vkr'' 

, 5 

m^ + (m + 2 )^ 

(m + 1 )ji 

W n2\ nL ) 


2(m + 1) 


(44) 


The minimum value of 


t is found from 
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which yields 


3/8, T \5/8/ \l/2 


f + \2/°, r 0 / 0 / , 

E W w 


(^ 5 ) 


This agrees with the data of Stein, Sanders, and Crate for ring- stiffened 
cylinders in torsion (ref. 8), at large values of Z^. It reduces to the 

isotropic solution for v = 0: 


| = 0.73l(t/R)^(-R/L)^ 2 


(k6) 


The result of equation (46) is 5-7 percent lower than the exact isotropic 
solution obtained by Batdorf for v = 0.313 • 


Long Cylinders 

For torsion on a long cylinder the classical solution for the iBO- 
tropic case is 


v = einfes + fj (W 

When this is substituted into equation (9) using equations (42), it is 
found that 



Since v = 0, then E/G = 2. Furthermore, 
(mitR/nL) « 1. Then equation (47) becomes 


it is assumed again that 
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2Tt _ (^f\(^) 2 

E ~ VmitR / \R/ 



Upon minimization of t with, respect to m and employing n 
the long cylinder, 


T 


1.754E 




1/4 


This reduces to the isotropic result, with v = 0, 

t = 0.272E(t/R) 3 / 2 


(49) 

= 2 for 

(50) 


(5D 
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TABLE 1 




THEORETICAL GENERAL INSTABILITY STRESSES FOR 
OKPHDTROPIC CIRCULAR CYILNEERS 


Loading 

Moderate-length cylinders 

Long cylinders 

Axial 

compression 

r 

a 2E 

(J a - 

c Et a 

P 2 ]^ + (GJ/2E) 

1/2 

r o i 1 / 2 

a 2E &%: + 

> 2 /t B ) + (E/2Gt) 

* c Rts (p 2 /t B ) + (E/2Gt) 

External radial 

or hydrostatic 

nrpHwnrv 

* 

Oy t= 5*511 

(o„\ J l Z / r .\ 

W IS; 


Oy - 3E(P f /R) 2 

Torsion. 

T . fcrf W 

\ “ / \ n ^y 

jL/2 

' 

* - 


a See equation (l4) for value of 0^. 
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